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摘要 : 氧气 在 质子 交换 膜 (PEM) 燃料 电池 阴极 的 还 原 反 应 特性 对 电池 性 能 有 重要 影响 。 本 文通 过 建立 电 催化 剂 Pt 表面 的 氧 
吸附 模型 , 及 氢 与 氧 在 Pt 表面 进行 的 氧化 还 原 反 应 模型 , 运用 分 子 动力 学 方法 模拟 研究 了 所 原子 与 氧 在 Pt 表面 的 反应 机 理 ， 分 
析 了 温度 对 氧气 还 原 反 应 特性 的 影响 。 研 究 发现 ， 氧 原子 和 和 氧 原子 在 Pt ge he 原 反应 的 控制 步骤 ， 随 着 
温度 升 高 ， 氧 气 还 原 反 应 的 速度 加 快 ， 但 温度 不 影响 氧气 还 原 反 应 各 个 步骤 中 的 产物 结构 。 研 究 结 果 对 加 强 理 解 PEM 燃料 电 
池 反 应 机 理 ， 推 动 燃料 电池 的 应 用 具有 重要 意义 。 
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Abstract: Oxygen plays an important role on Proton Exchange Membrane (PEM) fuel cell cathode reduction 
reaction. In this paper Oxygen adsorption model on catalyst Pt and the redox reaction model of hydrogen and oxygen 
on the surface of Pt are developed to study the reaction mechanism of hydrogen and oxygen in the Pt surface based 
on molecular dynamics simulation method. Analysis of the effect of temperature on the characteristics of oxygen 
reduction reaction is performed. The research results show that the first adsorption on the Pt surface for the hydrogen 
and oxygen atoms is the control step of the oxygen reduction reaction; with the increase of temperature, the oxygen 
reduction reaction becomes faster but the temperature does not affect the structure of the reduction product in each 
step. These results are very helpful in the understanding of PEM fuel cell reaction mechanism and promote the 
application of PEM fuel cell. 
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0 引言 两 步 四 电子 氧化 还 原 反应 过 程 站 。Adzic 等 中 认为 

百 : PEM 燃料 电池 阴极 氧气 还 原 反 应 以 四 电子 反应 机 理 
E 极 反应 尤其 是 阴极 氧气 的 还 原 反应 ， 对 PEM ”为 主 ， 二 电子 反应 机 理 与 四 电子 反应 机 理 同时 进行 。 
燃料 电池 的 性 能 有 重要 影响 。 近 年 来 ， 随 着 PEM 燃 ”Markovic 等 外 认为 二 电子 反应 机 理 更 适合 应 用 于 Pt 
料 电 池 技 术 的 发 展 ，PEM 燃料 电池 阴极 催化 层 氧化 和 Pt 合金 催化 剂 , 持 有 相同 观点 的 Sum 等 
还 原 反应 机 理 受 到 了 高 度 关注 。Pt 是 质子 交换 膜 燃 ”认为 二 电子 反应 机 理 比 四 电子 反应 机 理 更 适合 应 
料 电 池 目 前 最 常用 、 最 有 效 的 催化 剂 之 一 ， 因 此 发 ”用 于 氧化 还 原 反 应 。 根 据 反 应 过 程 中 在 能 量变 化 规 
生 在 Pt 面 上 的 氧化 还 原 反应 ， 成 了 研究 重点 之 一 。  ” 律 ，Liu 等 "基于 二 电子 理论 研究 了 氧化 还 原 反 应 机 
目前 ， 针 对 PEM 燃料 电池 阴极 中 氧气 和 质子 还 原 反 ， 理 ， 发 现在 第 一 步 反 应 时 释放 能 量 最 小 ， 在 第 四 步 
应 机 理 有 四 电子 反应 机 理 和 二 电子 反应 机 理 。 在 “反应 较 之 前 三 步 反 应 ， 放 出 最 多 的 能 量 。 在 第 一 个 
PEM 燃料 电池 中 ， 反 应 环境 是 弱酸 性 的 ， 当 Pt 中 杂 ”电子 传递 方面 ，Hartnig 等 ”WW 做 了 溶剂 分 
质 很 少 的 时 候 ， 氧 气 还 原 以 四 电子 反应 路 径 为 主 D。 子 的 重组 等 相关 的 研究 。 
在 弱 碱 性 环境 中 ，Miah 等 利用 伏 安 循环 技术 观测 到 
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Ana M. G6mez-Marin 等 外 认为 在 整个 反应 过 程 的 研 ee&6ss -ccd ps 
究 中 可 溶性 中 间 产物 的 减少 是 整个 氧化 还 原 反应 的 DD 

决策 步 。 催 化 剂 作用 下 ， 氧 气 和 质子 的 还 原 反应 得 
以 加 快 中 ,催化 剂 颗粒 的 大 小 和 分 布 ， 以 及 温度 和 酸 
碱 性 环境 也 对 其 反应 速度 和 机 理 有 重要 影响 nol nn 
为 了 加 快 PEM 燃料 电池 的 反应 速度 ， 各 种 催化 剂 及 


人 剂 都 被 研究 出 来 ， 并 进行 实验 对 比分 析 呈 1 “ 
[13 
966 
上 玉 文 献 主要 是 研究 PEM 燃料 电池 俊 化 剂 的 抽 a 
备 、 催 化 特性 、 电 极 特性 等 ， 尽 管 有 少数 文献 通过 68%e98os A 只 
模拟 分 析 了 燃料 电池 阴极 的 氧气 还 原 反应 机 理 ， 但 es sos 
基于 催化 剂 Pt 作用 下 的 氧气 还 原 机 理 以 及 影响 因素 ee 
的 研究 仍然 很 少 。 本 文 运用 分 子 动力 学 和 量子 力学 人 
的 基本 原理 ， 模拟 分 析 了 氧气 在 能 化 剂 pt 表面 还 原 Fig.1 O adsorption on the surface of the Pt 
= 有 反应 机 理 ， 分 析 了 温度 等 参数 对 Pt 催化 作用 下 氧气 
= ”还 原 反应 特性 的 影响 。 -20000 
S 1、 模型 -23836.5 四 
论文 模拟 采用 基于 密度 泛 函 理论 (DFT) 的 人 
CASTEP 进行 ， 分 别 模拟 高 温 PEM 燃料 电池 在 1 ] | 1 
一 393K、413K、433K、453K、473K 等 温度 条 件 下 ， zs 上 | ， 中 | 中 
本 ”氧气 背 压 为 0.1MPa 条 件 下 ， 发 生 在 阴极 催化 剂 Pt 和 
和- 表面 上 的 反应 ，Ox+4H_*2HO。 模 拟 中 的 结构 优化 -aaa 上 本 
全 ”采用 广义 梯度 近似 (GGA)， 交 换 关联 势 采 用 PBE 泛 
OQ 函 ， 计算 精 度 为 普 f 通 级 (Medium)， 计算 中 平面 波 斌 断 6 pn 克 10 区 3 本 的 
i 能 量 均 取 为 300eV， 自 浴场 循环 (SCF cycles) 设 为 Tmefs 
一 300， 精 ) 度 以 每 个 原子 能 量 收敛 至 2.0e-6ewatom 以 内 图 2 0; 吸附 于 Pt 表面 时 的 能 量变 化 
a 为 标准 。 分 子 动力 学 计算 采用 NPT 模式 ， 即 设 定 温 Fig.2 Energy with the O, adsorption on Pt surface 
ey 虑 、 时 日 页 ， ys 分 机 
= 反应 时 长 ， 固 定 晶 胞 质心 进行 分 子 动力 2、 反 应 路 径 
模拟 首先 对 Pt 原 胞 进行 结构 优化 ， 在 优化 后 的 在 弱酸 性 环境 下 ， 9 反应 路 径 为 ; 
Pt 原 胞 的 基础 上 切 出 Pt (111) 面 ， 构 建 4 层 Pt 原子 ， O,+4H' +4e 一 2H2O; 


添加 真空 层 厚度 为 2.0nm, 建立 1x2x1 的 超 唱 胞 。 约 二 人 
en et Pt 表面 的 氧气 还 原 反应 分 为 四 个 步 又 进行 。 第 
束 底 层 两 层 原子 ， 对 晶 胞 进行 弛 驳 ， 得 到 其 结构 参 个 SA A E 
Np 一 个 反应 步 又 是 将 第 仆 氧 原 卫 Hi 置 于 O01 与 0O;, 中 
数 为 : a=0.570905nm, b=0.988836nm, c=2.69921nm,， Dh 0 
a=B=Y=90”。 同 样 ， 对 0, 分 子 进行 结构 优化 ， 的 SS 2 


JE A *O+*OH; 《〈* 表 示 吸 附 于 Pt 表面 ) ， 如 图 3。 在 反 
测 得 优化 后 0 分 子 键 长 为 0.124014nm。 将 0; 置 于 ee 2 es 
Pt (111》 表面 的 真空 层 ， 0 将 吸附 于 Pt 表面， 对。 旋 过 各 中 ，H! 说 了 未 源 皖 汉 1 床 子 在 i 并 


吸附 的 结构 进行 优化 ， 得 到 能 量 最 低 的 0; 吸附 于 Pt Ne 

2 Lt 应 时 间 的 不 同 ， 可 形成 图 3(b)、(e) 所 示 两 种 生成 物 。 

CH1) 表面 的 稳定 结构 ， 如 图 1。 吸 附 过 程 中 能 量 A 训 于 po 化 的 

变化 如 图 2 所 示 。 能 时 随时 间 的 变 f 如 图 4 所 示 。 由 图 可 知 ， 在 反 

生成 反应 物 A 之 后 ， 会 有 一 个 瞬 态 能 量 稳定 期 ， 在 

此 0; 吸附 于 Pt (111) 表面 的 稳定 结构 作为 模 没有 得 到 充足 的 了 供应 的 情况 下 ， 最 终 分 解 为 中 间 

拟 计算 的 初始 模型 ， 运 用 分 子 动力 学 模拟 发 生 在 Pt 。 产物 B， 从 而 进一步 释放 能 量 。 但 实际 情况 是 了 的 

表面 的 氧气 还 原 反应 的 特性 。 计 算 所 得 到 的 Pt 晶 格 。“ 供 应 是 连续 而 充分 的 ， 因 此 第 一 个 步骤 中 的 反应 方 

常数 为 0.40nm， 该 常数 的 实验 值 为 0.392 nm， 在 误 程式 是 : *O0,+H 一 *OH, 而 不 会 出 现 *O0,H 一 *O+*OH 
差 范围 内 ， 该 计算 所 设 定 的 参数 是 合理 的 ， 所 得 数 。 ”的 反应 。 

据 在 有 效 范围 内 。 
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型 。 在 第 二 个 步骤 中 ， 加 入 第 二 个 氨 原 子 H， 发 生 

反应 : *O0,H+H 一 *O,Hs 一 *OH+*OH， 反 应 模型 如 图 
人 5。 反 应 中 将 H 原子 置 于 距离 01/、H1、0O, 等 同 距离 

0.3nm 处 ， 对 模型 进行 分 子 动力 学 模拟 。 可 以 发 现 : 

在 加 入 H; 之后， 了 逐渐 运动 趋同 于 0，,， 与 0; 成 
键 , 生成 x*O,H，， 并 最 终 分 解 为 两 个 *OH， 如 图 5(b)、 
图 5(c)。 


图 3 (a) 反应 物 
Fig.3 (a) the reactants C 


图 5(a) 反 应 物 
Fig.5 (a) the reactants 


3(b) 生 成 物 A 
Fig.3 (b) the product A 


图 5(b)H2 的 吸附 
Fig.5 (b) the adsorption of H2 


图 3(c) 生 成 物 B 
Fig.3 (c) the product B 
3 第 一 步 反 应 过 程 


《9? Fig.3 the first step reaction 


-23849 ee The reactant 

-23850 -RE 

ey 和 图 5(0O) 生 成 物 C 
-23852 ] Fig.5 (c) the product C 
23853 5 第 二 步 反 应 过 程 


Product B 


Fig.5 the second step reaction 


在 图 5(c) 即 生成 物 C 的 基础 上 ， 对 该 生成 物 进 


Energy /ev 


a] ， [PradactA| 行 优化 ， 并 加 入 第 三 个 氢 原 子 进行 第 三 步 反应 
0 100 200 了 | 500 600 700 800 2*OH+HH 一 *OH+H2O 。 初始 模型 如 图 6(a)， 将 H3 原 
Ime /Ss 


子 置 于 01 之 上 0.3nm 处 ， 之 后 对 该 模型 进行 分 子 动 
图 4 *O+H->*O2H->*O+*OH 体系 能 量 随 反应 时 间 变 化 力学 模拟 。0O1 在 与 Hi; 结合 后 ， 开 始 断 裂 与 Pt 原子 


Fig.4 Energy with the reaction of *O,+H—*O,H—*O+*OH 之 间 的 联系 ， 并 形成 水 4 分 子 脱离 离 ， 如 图 6(b)。 


日 于 实际 反应 时 H 是 连续 而 充分 地 供应 ， 故 将 最 后 对 第 三 步 得 到 的 产物 E 进行 结构 优化 ， ee 
生成 物 A 进行 结构 优化 ， 作 为 第 二 个 反应 步骤 的 模 ” ”入 第 四 个 氢 原 子 HL 作为 第 四 步 反 应 的 反应 物 ， 进 
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第 四 部 的 反应 : *sOH+H2O+H->2H2O。 第 四 步 的 初始 
反应 物 如 图 7(a) 所 示 ， 对 其 进行 分 子 动力 学 模拟 计 
算 ， 并 得 到 最 终生 成 物 F， 如 图 7(b)。 
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图 6(a) 反 应 物 
Fig.6 (a) the reactants 
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图 6(b) 生 成 物 E 
Fig.6 (b) the product E 


图 6 第 三 步 反 应 
Fig.6 the third step reaction 
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图 7(a) 反 应 物 
Fig.7 (a) the reactants 
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图 7(b) 生 成 物 F 
Fig.7 (b) the product F 
图 7 第 四 步 反 应 


Fig. the fourth step reaction 
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模拟 设 定 压强 等 同 、Pt 表面 状态 相同 和 放置 了 
原子 位 置 相同 , 分 别 分 析 温 度 在 393K、413K、433K、 
453K 和 473K 时 对 反应 中 四 个 步骤 的 影响 。 反 应 的 
四 个 步骤 是 : 


*O,+H—*O,H GD 
*O2H+H-~*O2H2 一 *OH+Y*OH © 
2*O0H+H*OH+H,O G) 
*OH+H2O+H-~2H2O 由 


每 个 反应 步骤 都 是 在 前 一 个 步骤 得 到 生成 物 的 基础 
上 ， 在 相同 的 位 置 加 入 反应 物 H 原子 进行 下 一 步 的 
反应 。 


图 8 是 当 温 度 条 件 为 433K 时 , 四 个 步 又 反应 过 
程 中 的 能 量 随时 间 的 变化 。 在 HH 的 初始 位 置 距离 氧 
原子 等 同 的 情况 下 ， 每 一 步 反 应 所 需 的 时 间 不 同 ， 
We 一 步 反 应 所 占 时 间 最 长 ， 达 到 130fs 左右 。 第 

、 四 步 反 应 所 需 时 间 均 小 于 此 值 ， 三 者 所 需 
反应 时 间 均 小 于 50fs。 由 此 可 见 ， 在 Pt (111) 表面 


进行 的 O, 与 H 的 反应 中 ,HH 在 0， 上 的 首次 吸附 所 
需 的 时 间 在 整个 反应 过 程 中 是 最 长 的 ， 这 意味 着 这 
个 步骤 是 整个 反应 中 的 控制 步骤 。 


Energe/eV 
Dn 
3 
Oo 
名 
[= 


0 50 100 150 200 
Time/fs 


Fig.8 four steps in the process of system energy change over time 


fewer than 433 Kelvin temperature 


图 9 是 整个 反应 控制 步骤 第 一 步 反 应 速度 随 温 
度 的 变化 曲线 。 图 中 显示 : 随 着 温度 的 提升， 反应 
成 键 之 前 所 需 的 时 间 逐 渐 减 少 ， 即 分 子 运 动 的 速度 
随 温度 的 升 高 而 增加 。 同 时 ， 随 着 反应 温度 的 提升 ， 
成 键 后 达到 稳 态 的 时 间 也 显著 缩短 ， 即 反应 的 速率 
随 温 度 的 升 高 而 加 快 。 

表 1 统计 了 各 温度 条 件 下 控制 步 又 第 一 步 反 应 
O 和 Hi 吸附 瞬间 H 原子 的 瞬间 速度 以 及 吸附 前 卫 
原子 的 平均 速度 , 0; 和 Hi; 之 间 成 键 的 时 间 和 稳定 后 
结构 的 最 优 键 长 。 


Energe/eV 


图 9 不 同 温度 下 第 一 步 反 应 的 时 间 和 能 量 的 变化 


different reaction temperature 


表 1 不 同 温度 下 的 Pt 原子 间距 及 0, 与 Hi 的 成 键 时 间 


Table 1 the distance between two Pt atoms and the bonding time 


between O2 and H1 for different temperature 


温度 393K 413K 433K 453K 473K 
吸附 前 H 平 ”0.000272 0.000374 0.000731 0.001054 ” 0.001111 
均 速 度 08 595 588 976 557 
/0.1nm-fs’! 
成 键 时 间 /fs 229 186 113 gb 66 
O: 与 HI 最 优 0.990 
0.990 0.988 0.988 0.991 
键 长 /0.1nm 9 9 


综合 表 中 数据 可 以 发 现 ， 随 着 温度 的 上 升 ，H 


原子 吸附 前 和 吸附 瞬间 的 运动 速度 显著 增加 ， 成 键 
需要 的 时 间 缩 短 ， 但 温度 对 生成 物 稳定 状态 的 最 优 
键 长 影响 很 少 ， 可 以 忽略 不 计 。 
4、 结 论 

本 文 根 据 密 度 泛 函 理论 ， 对 Pt (111) 表面 的 
0; 的 还 原 反 应 进行 分 子 动力 学 模拟 计算 ， 根 据 计算 


的 结 


果 ， 可 得 到 如 下 的 结论 : 


1、 在 Pt(111) 表面 上 进行 的 O, 与 再 的 反应 过 
程 中 , 如 果 了 原子 能 持续 和 充分 供给 ,其 反 
应 步 又 中 的 第 一 步 反 应 只 能 是 : *O0,+H 一 
*O,H。 


2、 在 Pt (111) 表面 上 进行 的 0, 与 H 的 整个 反 
应 过 程 中 ， 第 一 个 步骤 ， 即 Hi 与 0; 的 吸附 
时 间 最 长 ， 是 整个 反应 的 控制 步骤 。 

3、 随 着 温度 的 升 高 ，H 原子 的 运动 速度 加 快 ， 
和 0 吸附 成 键 时 间 缩 短 ， 表 明 温 度 升 高 可 
以 提高 反应 速率 , 但 不 会 影响 各 个 反应 步骤 
中 生成 物 的 结构 。 
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